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ABSTRACT 

We investigate the formation of GMCs in spiral galaxies through both agglomeration 
of clouds in the spiral arms, and self gravity. The simulations presented include two- 
fluid models, which contain both cold and warm gas, although there is no heating or 
cooling between them. We flnd agglomeration is predominant when both the warm and 
cold components of the ISM are effectively stable to gravitational instabilities. In this 
case, the spacing (and consequently mass) of clouds and spurs along the spiral arms is 
determined by the orbits of the gas particles and correlates with their epicyclic radii (or 
equivalently spiral shock strength). Notably GMCs formed primarily by agglomeration 
tend to be unbound associations of many smaller clouds, which disperse upon leaving 
the spiral arms. These GMCs are likely to be more massive in galaxies with stronger 
spiral shocks or higher surface densities. GMCs formed by agglomeration are also found 
to exhibit both prograde and retrograde rotation, a consequence of the dumpiness of 
the gas. At higher surface densities, self gravity becomes more important in arranging 
both the warm and cold gas into clouds and spurs, and determining the properties of 
the most massive GMCs. These massive GMCs can be distinguished by their higher 
angular momentum, exhibit prograde rotation and are more bound. For a 20 M© 
pc~^ disc, the spacing between the GMCs flts both the agglomeration and self gravity 
scenarios, as the maximum unstable wavelength of gravitational perturbations in the 
warm gas is similar to the spacing found when GMCs form solely by agglomeration. 
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1 INTRODUCTION 

The accumulation of the ISM into giant molecular clouds 
(GMCs) represents the earliest stage in star formation. The 
properties of GMCs, such as turbulence and magnetic fleld 
stre ngth regulate how star formation evolves o n local scales 
(e.g. lMcKeellT999l : IPudritd [200i : iLarsonl liooi ) . and are in- 
trinsically linked to many of the issues in star formation, 
such as the time for stellar collapse (e.g. (Elmegreen 2007)). 
Thus understanding how GMCs form, and how their char- 
acteristics depend on the dynamics of galaxies and nature 
of the interstellar medium (ISM), is essential for progress in 
star formation. The formation and evolution of GMCs is also 
interrelated to the global properties of the ISM. For exam- 
ple, GMC formation by coalescence is much more likely in 
a cold, clu mpy medium, con i pared to a warm diffuse en- 
vironment (|Dobbs Bonneli l2006h , where instabilities in 
the ISM are required. In turn stellar feedback controls the 
ejection of hot gas back into the ISM and generates tur- 
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bulence fe .g. iMac L ow & Klessen 20^; Ijoung Mac Lowl 
2006; de A villez k Breitschwerdt 200'$^. 



There have been numerous suggest ions of how GMCs 
form (as described in a recent review bv lMcKee Ostr ikerl 
I2OO7I and references therein). However they predominantly 
fall into two categories: either GMCs form by the agglom- 



the 



eration of smaller clouds of g as (e.g. Field Sasla\yl 

1965; TaflF & Savedog Il972l : IScoville HershI Il979l : 
ICasoli CombesI 11983 ). or through instabilities 
ISM. The latter incl u de gr av itational (ICowid 
Elmegreen ElmegreenI Il983l : iBalbus Cowid 
La Vigne et al.1 l2006l\ magnetic (either Parker msta- 



1981 



1985 



bilities (iMouschovias eraD 1 1 9741 : iBlitz Shul ll98Qh or 

MRI (|Kim et al.1 l2003h) or therm al instabilitie s ~"(lFieldl 

Il965l : iKovama Inutsukal I2OOOI : IStiele et al.1 I2OO6I I. 

Another recent suggestion is that GMCs form in col- 

I 1 1 \ 

liding or turbulen t flows (|Vazquez- Semadeni et al. 1995; 

Balle s teros-Paredes e ^^ Il999l : ^Izquez-Semadeni et al.| 

I2OO6I : iHeitsch et al.ll2006l ). This however requires a mech- 

anism to produce these flows, which is most likely gravity 



Elmegreen et al.l 2003ll. spiral shock s 



Kovama &; Inutsukall200"ol : iBergin et"alll2004l ) 



or super novae 
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Coalescence of smaller clouds was first instigated to ex- 
plain GMC form ation., but the timescales of 10^ Myr for 
formation (Kwan 1979) were presumed too long; to accoun t 
for the observed properties of GMCs telitz Shul Il98(t >. 
Furthermore star formation may be expected to disrupt the 
const ituent clouds before a m ore massive GMC is assem- 
bled (jMcKee k Ostrikedl2QQ7l ). By including spiral density 
waves, the ti me f or formation by coalescence is reduced to 
IP's of Myrs dCasoli k Combes.. 1982. : .Kwan k Valdes..l987 : 
[Roberts k Stewartlll98"7E 



In more recent reviews, gravitational instabilities gen- 
erally a£2ea£_the preferred m echanism for GMC formation 
(e.g. lElmegreenI 1990: McKe e k Ostrikerl l2QQ7l ) . Theoreti- 
cal analysis (lElmegreenI 1198^ 1 and numerical simulations 
(|Kim et al.1 l2QQlL I2002l 7 suggest that the Parker instabil- 
ity alone produces insufficient density enhancements, and 
GMCs can form by gravitational instabilities in relatively 
short timescales (Elmegreen 1990). Predictions from gravi- 
tational analysis have also been compared to the observed 
masses and sizes of GMCs. In spiral galaxies, GMCs ap- 
pear to be regularly spaced along the spiral arms. This 
spacing, and the mass of GMCs can be interpreted in 
terms of gravitational perturbations to the gas (Cowie 
19811: lElmegreen k Elmegreenlll983l : iBalbus k Cowielll985 



La Vigne et al.ll2006l ). For a sound speed of 7 km s ^, the 



estimated mass of a GMC is 10^ - 10^ M©, whilst the 
characteristic spacing is expected to be around 3 times the 
width of the s piral arm, c orresp onding approximately with 
observations (jElmegreenl ll995l V Finally, GMC formation 
by gravitational instabilities h as further been endor s ed by 
recent numerical s imulations (|Kim k Ostrikerl |2002| , l2006l : 
IShettv k Ostrike r 2006) which produce spacings and masses 
in line with the theoretical predictions. 

The problem of molecular cloud formation has 
been reopened in the last decade, with computational 
resour ces now enabling nu n ierical simulations on bo th 
local (iKim k Ostrikeil I2OO2I: iGlover k Mac Lowl l2007a 



IVazauez-Semadeni et al.1 '2006*: 'Heitsch et al. 2006^ and 
galactic ([Shettv k Ostriker 2006; Dobbs et al. I2OO6I : 
iDobbs k Bonnelll I2OO8I : ["asker et all I2OO8I ) scales. These 
simulations show the fo r matio n of GMCs by self g ravity 
(iKim k Ostriked I200I I2OO6I : IShettv k Ostriked I2OO6I : 
iGlo ver k Mac L owl l2007al l, turbulent or colliding flows 

(IVazauez-Semadeni et al. r_^Qo&_ ■Glover k Mac Low 



2007bl: iHeitsch et al.1 I2OO6I : iHennebelle et al.1 1200^ 
Heitsch et al. I "2008"), combined Parker and thermal 



stabilities (jKosihski k Hanasz 2007^) , and Kelv in Helmholtz 
instabilities (Wada k Koda 2004; Wada 20 03) . In pre vious 
results (Dobbs k Bonnell 2006 : Dobbs et aD I2OO6I ). we 
showed that molecular clouds, and inter-arm spurs, can 
form as a result of cold gas passing through a spiral shock. 
The formation of this structure occurs without self gravity, 
and is still present when magnetic fields are included, 
although magneti c pressure acts to sm ooth out clumps in 
the spiral shock (jDobbs k Pried I2OO8I ) . The formation of 
clouds in these calculations bears most resemblance to the 
collisional models. Clumps of gas are forced together by 
the spiral shock, and agglomerate into larger structures, 
which become most discernible as interarm spurs extending 
from the spiral arms. Generally in simulations of grand 
design spirals, the spiral pattern is assumed to be long-lived 



(incl uding the present work), although this ni ay not be the 
case (jMerrifield et al.ll2006l : fshettv et al.ll2007l ). 

This paper concentrates on two possibilities, the forma- 
tion of GMCs by agglomeration, and by gravitational insta- 
bilities. In particular, we discuss the relative contribution 
to GMC formation for different initial conditions. When the 
surface density is lower, GMCs form by the agglomeration of 
smaller gas clouds as their orbits converge in spiral shocks, 
similar to the model by Roberts k Stewar t ( 1987 ). In this 
scenario, the spacing and size of clumps depends on the time 
(or equivalent ly distance) gas spends in the spiral arm, and 
interactions between clumps during orbit crossings. We show 
also the formation of GMCs form by gravitational instabil- 
ities for a high surface density disc with solely warm gas, 
and show that for a high surface density disc containing 
cold gas, both agglomeration and gravitational instabilities 
contribute to GMC formation. We assess the differences in 
structure of the disc, in particular the separation of spurs, by 
performing MHD calculations with and without self gravity, 
varying the strength of the potential, disc mass, temperature 
and magnetic field strength. 



2 CALCULATIONS 

We use the 3D SPMHD c ode, a ver sion of the S PH code 
originally written by Benz (|Benz et al . 1990) and B ate et aP 
(|l995l ) . The code has been extended to include magnetic 
fields (jPrice k MonaghanI liooHI . l2004l ). For the simulations 
presented here, the magnetic field is represented by Eu- 
ler potentials (see also Section 2.2). The code also uses 
a variable smoothing length, such that the density p and 
smoothing length h are solved iteratively according to 
IPricefcMonaghanldioOTi ). Artificial viscosity is included to 
treat shocks, wit h the standard values a = 1 and (3 = 2 
(|Mona ghan 1997). 

Long range gravit ational forces are calculated using a 
binary tree algorithm (|Benz et al.lll99(3h . Over short scales, 
the gravitational forces are softened by setting the soften- 
ing length equal to the s moothing length, as described in 
'Price Mona ghanI (|200 7'l. This method conserves both mo- 
mentum and energy. Tests of the formulation of gravitational 
forces, and the implementation of gra vitational softening are 
included in lPrice k Mon aghan (^20071). The code and the im- 
ple mentation of Euler p otentials are describ ed in more detail 
in lPrice k Batd (|2007l i and lDobbs k Pried (|2008l ). 



2.1 Thermal distribution of the gas 

The calculations presented here comprise of single phase and 
'two-fluid' models. The gas in the single phase calculations 
is either cold (100 K) or warm (10^ K) whilst in the two 
fluid model, half the particles are allocated a temperature 
of 100 K and half 10^ K. In aU cases, the calculations are 
isothermal, thus although the two fluid model includes cold 
and warm gas, there is no phase transition between the two 
components. 

This two fluid model is clearly a very simplified ap- 
proach with numerous limitations. In particular, this model 
assumes that cold gas is widespread in the ISM, and enters 
the spiral shock. Recent calculations indic ate that gas tend s 
to heat up after passing a spiral shock ( Kim et al] I2OO8I I , 
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in which ca se the gas e ntering the shock may be warm. 
However in Wadal 120081), there appears to be a consider- 
able degree of cold gas in the midplane of the disc (Fig. 2) 
whilst some observations and theoretical analy sis suggests 
cold HI or H2 may be entering spiral sho cks (|Vogel et al.l 
ll988l : IPringle et al.ll200l] : lGibson et al.ll2007l l. I n calculations 
which include thermodynamics of the ISM (jPobbs et al.l 
12008'^. we find that although gas entering the shock is usually 
warm, it tends to cool very quickly and produces clumpy, 
cold medium where agglomeration can occur. Future sim- 
ulations with full MHD will be able to address this issue 
further. 

In the two fluid models, the warm and cold gas gener- 
ally separate out during the simulations, as the spiral shock 
induces much higher densities in the cold gas. In the agglom- 
eration scenario, these results are not dependent on using a 
two fluid medium, since the structure in the cold gas is es- 
sentially the same as when only cold gas is used. In the low 
surface density case, where there is little substructure in the 
warm gas, and the cold clumps are largely separate from 
the warm gas, the pressure from the warm phase confines 
clumps and spurs to higher densities. Hence these features 
are easier to distinguish in analysis of the two fluid models. 
It is also useful to compare the structure of the disc for a 
given surface density when the gas is cold, warm and a mix- 
ture of both, and the latter is the most realistic distribution 
for the ISM of the 3 we investigate. 



2.2 Initial conditions and details of simulations 

We model an isothermal galactic disc between radii of 5 
and 10 kpc. The gas is initially distributed uniformly with 
z < 150 pc for the single phase cold calculations, and 
z < 400 pc for the two fluid and warm calculations. With 
time, the scale height is t ypically 20-100 pc fo r the cold gas 
and 300 pc for the warm (jPobbs Pricell2008l ). We perform 
calculations with 3 different surface densities (see Table 1), 
and for the majority of calculations, we use 4 million parti- 
cles. The corresponding mass of the discs are 10^ 2 X 10^ 
and 4 x 10^ M©. This gives a mass resolution of 250, 500 
and 1250 M© for the 4, 8 and 20 M© pc~^ surface den- 
sity calculations respectively. In the appendix, a resolution 
study is included, with simulations of 1 and 8 million parti- 
cles. All particles have the same mass, hence the two fluid 
simulations contain the same mass of warm and cold gas. 

The velocities in the plane of the disc follow a rotation 
curve corresponding to the disc component of the potential 



described in terms of Euler potentials by 



1 2, 

2^0 log 



2 2 
r ^ z 

Zq 



+ 1 



(1) 



where Rc=l kpc, vo = 220 km s~^, and Zq=0.7 kpc is a mea- 
sure of the disc scale height. This produces an essentially flat 
rotation curve for the radii over which the particles are dis- 
tributed. We also impose a velocity dispersion by selecting 
velocity perturbations from a Gaussian of mean 6 km s~^. 
The velocity dispersion constitutes the z component of the 
velocities for the particles. 

The magnetic fleld in these simulations is initially 
toroidal and of uniform strength. A toroidal fleld can be 



I3e 



where r 



1 2 



(2) 
(3) 

+ + z'^ and — cos~^ {z / r) . The relative 
strength of the magnetic fleld is given by the mean plasma 
beta, shown in Table 1. This is deflned as the ratio of gas to 
magnetic pressure 



2AtopoC 
IBP 



(4) 



where po is the average density of the disc. We stress that 
although the magnetic fleld is uniform initially, there are 
large variations locally in the fleld strength and /3 as the 
simulation progresses. 

During the simulation, the gas is subject to an ex- 
ternal potential. The pot ential includes a galactic halo 
(jCaldwell Ostrikerl ll98lT ) and the logarithmic disc com- 
ponent given in Eqn 1. T he potential also incorporates a 4 
armed spiral perturbation ( Cox G6mezll2002l : lDobbs et al] 
.2006 ). with a pattern speed of 2 x 10^ rad yr~^, and a 
pitch angle of 15"^. Calculations are performed with different 
strength potentials, as shown in Table 1. The strength F is 
determined by the maximum radial force of the perturba- 
tion compared to the underlying disc potential, at a radius 
of 7.5 kpc, i.e. 



dr 



' di/jdi. 



dr 



(5) 



whe re V^gp is the p otent ial for the spiral perturb ation (see 
also'Roberlsl (|l969h and lShettv Ostriked (|2006l ll. We run 
the simulations for at least 250 Myr, or until runaway col- 
lapse occurs. 



2.3 Calculation of Q, Xmax and ve 

This section describes the calculation of the parameters Q, 
Xmax and rE listed in Table 1. These parameters are dis- 
cussed in relation to the origins of GMCs in our simulations 
in Section 3.1.2. As mentioned in the introduction, the sepa- 
rations and masses of GMCs in galaxies can be compared to 
those predicted by theoretical models. We compute the ex- 
pected separation and mass of large scale structures result- 
ing from gravitational instabilities, as well as the st ability of 
the d isc, as measured by the Toomre parameter Q (jToomrd 
1 19641 ). However these are only approxima tions as we d o not 
include a correction for the spiral shock ( Balbuslll9 88V and 
use the total column density, even in the two fluid models. 
We discuss any differences from these assumptions below. 
The Q parameter is calculated by 



Q — kCs/ttGH 



(6) 



where n = V2vc/R {vc = 220 km s""*^ is the rotational veloc- 
ity and the radius assumed is R=7.5 kpc), and — ksT / ji 
where /x = 1.4 is the mean molecular weight in these simu- 
lations. The disc is considered stable if Q > 1. For the two 
fluid models, we take the total surface density, but use the 
sound speed of the cold and warm gas to provide parame- 
ters Qc and Qu (Table 1). This is probably reasonable for 
the high surface density case where structures containing 
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Two phase calculations 



Model T (cold) E (3 Self F Qc Spacing r^; Xmaxc >^maxh ^dump 

(K) (M0PC-2) gravity (%) (kpc) (kpc) (kpc) (kpc) (10^ M©) 



A 100 4 4 Yes 2 0.5 5 0.84 ±0.04 0.45 0.066 6.6 0.5 

B 100 4 4 Yes 4 0.5 5 1.05 ±0.06 0.53 0.066 6.6 1 

C 100 4 0.4 Yes 4 0.5 5 0.91 ±0.04 0.53 0.066 6.6 0.2 

D 100 4 4 Yes 8 0.5 5 1.52 ±0.09 0.73 0.066 6.6 3 

E 1000 4 4 Yes 8 0.5 5 1.33 ±0.09 0.73 0.21 6.6 0.7 

F 100 4 4 Yes 16 0.5 5 1.66 ±0.04 0.85 0.066 6.6 4 

G 100 4 4 No 4 - - 1.04 ±0.04 0.53 - - 0.4 

H 100 8 4 Yes 4 0.25 2.5 1.12 ±0.07 0.53 0.033 3.3 3 

I 100 20 4 Yes 4 0.1 1 1.26 ±0.06 0.53 0.013 1.3 20 



Single phase calculations 



Model 


T 


E 




Self 


F 


Qc 




Spacing 




^maxc 








(K) 


(M0PC-2) 




gravity 


(%) 






(kpc) 


(kpc) 


(kpc) 


(kpc) 


(106 Mo) 


J 


100 


4 


1 


Yes 


4 


0.5 




1.00 ±0.03 


0.53 


0.066 




1 


K 


10000 


4 


4 


Yes 


4 




5 




0.53 




6.6 




L 


100 


20 


4 


Yes 


4 


0.1 






0.53 


0.013 






M 


10000 


20 


4 


Yes 


4 




1 


1.43 ±0.06 


0.53 




1.3 


8.5 


N 


10000 


32 


4 


Yes 


4 




0.63 


1.04 ±0.04 


0.53 




0.8 
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Table 1. Table showing the parameters used in these calculations. All calculations use 4 million particles and in the two phase calculations, 
half the gas is distributed in the temperature indicated in the table, and half is 10"^ K. (3 is the ratio of thermal to magnetic pressure 
for the cold gas, and F is a measure of the strength of the spiral perturbation, r^ is the epicyclic radius at R = 7.5 kpc and Q is the 
Toomre parameter. The calculation of these quantities is described in Sections 2.2 and 2.3. The spacing is calculated (also described in 
the text) after 265 Myr for all models except I (with time 130 Myr), L (time = 180 Myr) and N (time = 85 Myr) where gravitational 
collapse prevents the simulations continuing. The final column is the mass of the most massive clump in a 5 kpc x 5 kpc section of the 
disc, using a clump-finding algorithm with a surface density threshold of 5E. In models K and L, regular large scale GMCs and spurs do 
not form within the time of the simulation. 



both cold and warm gas form along the spiral arms. How- 
ever this treatment probably underestimates Q in the low 
surface density cases (where the fluids are more separate) 
and the disc is more stable than suggested in Table 1. This 
calculation of Q also neglects magnet ic fie lds, which would 
modify Q by a factor of 1 + 1//3 (e.g. [Sh3l^92). 

For a thin disc, the dispersion relation arising from sta- 
bility analysis of a differentially rotating disc is 



- /^^ + 2iTGi:\k\ - kci = (7) 

fe.g. lBinnev &: Tremaine|[l987l V The separation of the clouds 
is expected to be approximately 



= 2cJ/GE, 



(8) 



the wavelength corresp onding to the peak g rowth - 
rate of pe rturbations (jElm egreen Elmegr eenI llQSSi : 
iBalbus Cowie 1985^) whilst the minimum length at which 
perturbations become unstable for a thin disc, the Jeans 
length, is c^/GS. The mass enclosed within a cloud is 



M 



(9) 



the Jeans mass for a thin rotating disc. Again S includes 
both the warm and cold gas in the two fluid models. How- 
ever if only one component is included, this still produces 
a very different spacing from that of the low surface mod- 
els. Also if we take the surface density of the spiral arms, 
as oppose to the average value, Xmax decreases still further 
for the cold gas. The warm gas does not experience such a 



strong increase in density in the spiral arms, but our \max 
may be a high estimate. Overall though, Xmax turns out 
to be very different from the spacing in calculations where 
gravitational instabilities are not believed to be dominat- 
ing the large scale structure, but similar when gravitational 
instabilities are thought to be responsible. 

For calculations where the spacing is believe to corre- 
spond to formation by agglomeration, we compare the spac- 
ing to the epicyclic radius of the stellar orbits, a measure of 
the strength of the spiral potential. The epicyclic radius is 
calculated for a 2D test particle simulation subject to the 
spiral potential alone, at a radius of 7.5 kpc. Particles are 
placed in a disc, and those with an average radius between 
7.4 and 7.6 kpc over 2 rotational periods are selected. The 
epicyclic radius is then the average (r max T^min ) /2 of these 
particles (although this is a radial measurement in the disc 
as oppose to along a spiral arm) . 



3 RESULTS 

We show the structure of the galactic disc for several simu- 
lations in Figs 1, 2 and 7. Except for models K, N (no cold 
gas) and G (no self gravity), both agglomeration and accre- 
tion by gravity contribute to the formation of GMCs. When 
S ^ 8 Mq pc~^, clouds and spurs form predominantly by 
agglomeration. However in models I, M and N, when S = 20 
M0 pc~^, the gas becomes gravitationally unstable, and self 
gravity has a greater effect on the structure of the disc and 
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the properties of GMCs. Thus we discuss the different sur- 
face density regimes separately in the next sections. 

3.1 Formation of clouds by agglomeration (low 
surface density models) 

In Fig. 1, the galactic disc is displayed for models B, G, and 
D, which compare results with and without self gravity (a 
and b), and with a stronger potential (c). All three assume 
a two fluid medium. The last panel (d) shows the disc in 
model J, where there is only cold gas. We see from b) and 
d) that the structure in the disc is dominated by the cold gas, 
rather than the warm, with essentially the same structure in 
the cold gas and two fluid simulations (J and B). The warm 
gas in these models does not experience significant density 
perturbations along the spiral arms. Fig. 2 shows a disc with 
the same surface density, but all the gas is warm (run K). 
In this case, the disc is stable to gravitational instabilities, 
and the gas too warm for clump agglomeration to occur, so 
there is no substructure. 

The top two panels of Fig. 1 compare simulations with 
(a) and without (b) self gravity. When there is no self gravity, 
spurs are still evident. With self gravity, the gas is arranged 
into more coherent clumps and spurs. However, the large 
scale structure, in particular the largest spurs which extend 
between the spiral arms, is very similar in both cases. Look- 
ing at a more detailed section of spiral arm (Fig. 3), there 
is even a 1 to 1 mapping between the spurs, with and with- 
out self gravity, although there is more fragmentation of the 
spurs without self gravity. 

We interpret the formation of clumps and spurs in these 
simulations as due to agglomeration of gas in the spiral 
shock, rather than self gravity. In Fig. 4, we show the gas 
which constitutes one of the clumps in model B. The clump 
in the lower panel forms by the agglomeration of smaller 
clumps during the spiral shock passage. The same process 
occurs regardless of whether we include self gravity. How- 
ever it is easier to distinguish spurs and spiral arms clumps, 
when self gravity is included. 

By viewing a movie of gas flowing through a spiral 
arm, it is evident that clump agglomeration is not partic- 
ularly efficient. Clumps do not always collide since the size 
scale of the clump is less than the width of a spiral arm, 
and collisions result in fragmentation as well as agglomer- 
ation. Se lf gravity increases th e chances of agglomeration 
(see also iKwan k Valdeslll987l ) and leads to more distinct 
clumps (Figs 1, 3 and 4). It is apparent from Fig. 4 that the 
time to assemble a GMC structure from constituent clumps 
is around 50 Myr, i.e. comparable to the time for a spiral 
arm passage. 

We describe the process of agglomeration in more de- 
tail, as well as explaining why the spacing between clouds 
and spurs agrees with agglomeration as the main mecha- 
nism behind their formation in Section 3.1.2. Two further 
issues with the agglomeration model, the nature of the gas 
prior to the shock, and the dependence of cloud properties 
on numerical resolution are discussed in the appendices. 

3.1.1 Wiggle instability 

An alternative explanation, which also does not require 
self gravity or magnetic fields, is that clumps and spurs 



are formed by the wiggle ( Kelvi n- Helmholtz) instability 
(|Wada k Kodal l20Q4l : IWadal l2008l ). Although this a pos- 
sibility, we favour the agglomeration scenario as respon- 
sible for the large scale structure in our models. We 
still obtain structures at low resolution even though the 
Rey nolds numbers are very low for resolving K-H instabili- 
ties (iDobbs k Bonnellll2006l ). At lower resolution the insta- 
bility appears weaker (fig. 4, Wada k Koda (2004)) whereas 
spurs and clumps are actually more distinct in our sim- 
ulations (see Appendix B). Secondly, the rotation curve 
of our disc is fiat over the radii modeled, hence the disc 
should be rela t ively stable to K-H instabilities (although 
IWada k Kodal (|200 4) point out that spiral streaming mo- 
tions alone may induce instability). Thirdly agglomeration 
and merging of clumps is evident in simulations - even if KH 
instabilities are present initially, the gas evolves into clumps 
and the dynamics are then dominated by interactions be- 
tween the clumps. Finally, as will be described in Section 
3.1.2, we find a correlation between the spacing of clumps 
and the strength of the spiral potential, which is consistent 
with formation by agglomeration. It is not clear how the 
spacing varies if the structure arises from the wiggle insta- 
bility. 

3.1.2 Separation and masses of clouds 

In this section, we argue that the separation of clouds in 
the low surface density models is consistent with formation 
by agglomeration rather than gravitational instabilities. Ta- 
ble 1 lists the separations between the main clouds/spurs 
along the spiral arms. These separations are calculated by 
selecting particles from a section of spiral arm. The particles 
are then transformed to co-ordinates parallel and perpen- 
dicular to the spiral shock, and particles binned into 5 pc 
sections over the coordinate parallel to the shock (thus the 
spacing includes the curvature of the spiral arm). We used 
two methods to obtain the spacing. Firstly, the average dis- 
tance between the peaks along the spiral arm was calculated 
by selecting the main peaks by eye. This was repeated for 
the same region over each of the 4 spiral arms, and the av- 
erage from the 4 spiral arms listed in Table 1. The error is 
also calculated from the spread in these values. We also ap- 
plied Fourier analysis to the distribution of particles along 
the arm. For the models with stronger potentials (F=6% 
and 12%), the average spacing over the 4 arms is similar to 
those listed in Table 1, although with more spread. However 
for the lower strength potentials, the Fourier transform for 
a section of spiral arm sometimes produced peaks of similar 
magnitude. For these results it was inappropriate to extract 
a value for the spacing. Nevertheless we show in Fig. 5 some 
sample Fourier transforms where a reasonable signal was re- 
covered for different strength potentials. 

For the low surface density models, the separations pre- 
dicted from Jeans analysis of the warm gas are much too 
large compared to those in the simulations. The difference is 
not surprising since the warm component is gravitationally 
stable due to the low surface density. However, as evident 
in Table 1, the cold gas is expected to collapse. Due to the 
dependence of the length scale in Eqn (8) on sound speed, 
the size scale of instabilities is only lO's of parsecs in the 
cold gas, and is therefore much smaller than the separation 
between the largest clumps. Thus self gravity is important 
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Figure 1. The galactic disc is shown for models B (a), G (b), D (c) and J (d) after 265 Myr. The top 2 panels compare results with (a) 
and without (b) self gravity. The spurs and structures in the arms are clearer when self gravity is included, although the overall structure 
is not significantly different. Both the lower panels (c and d) include self gravity. For (c), a stronger potential is used, and the spurs and 
clumps in the spiral arms are spaced further apart. The fourth panel (d) is from model J where all the gas is lOOK. In all cases the disc 
surface density is 4 pc~^. 



in the cold gas, but does not determine the large structure. 
The smaller spherical clumps shown in Fig 3, of size scales 
50 pc or so, are comparable to the size scale associated with 
gravitational instabilities, and it is on these scales self grav- 
ity is having a more significant effect, and produces more 
coherent structures. The stability of the cold gas in the disc 
is discussed further in Section 3.1.3, and in particular why 
the cold gas is effectively stable and global collapse does not 
occur. 

The spur separation is however found to increase with 
the strength of the spiral perturbation, or equivalently 
strength of the spiral shock (Table 1 and Fig. 5). This be- 
haviour is contradictory with the Jeans analysis. If the spi- 
ral perturbation increases, we would expect Xmax to stay 
the same or decrease, depending on whether we take the 
average surface density, or the sur face density along the spi- 
ral arm. IShettv &: Ostriked (|2006l ) find that the spacmg re- 



sulting from gravitational instabilities is similar for different 
strength potentials. We also find that the spacing decreases 
for a calculation with 1000 K gas (E), corresponding to a 
decreasing shock strength, compared to that where the cold 
component is 100 K (D). If cloud formation is due to grav- 
itational instabilities, the spacing would be expected to in- 
crease due to dependence on (Eqn. 8). Finally, when we 
double the disc mass, the spacing does not change signifi- 
cantly, as would be expected if the spacing is determined by 
the orbits. On the other hand if self gravity is dominant, we 
would expect the spacing to halve. 

Rather than self gravity determining the size and spac- 
ing of the spurs and clouds, we relate these properties to the 
convergence of orbits through the spiral shock. With grav- 
itational instabilities, when a density perturbation occurs, 
there is a surrounding region (governed by the wavelength 
of the perturbation) of gas which will be subject to col- 
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Figure 2. The column density of the disc is shown for model 
K after 235 Myr. All the gas is warm and the surface density 
is 4 Mq pc~^. There is no substructure in the disc since the 
warm component is stable to gravitational instabilities, and the 
gas pressure prevents any clumpy structure evolving. 

lapse. Analogous to this, in our models there is a locus of 
gas particles whose orbits will pass close to the initial per- 
turbation. These gas particles are on slightly different orbits, 
so will have different velocities and angular momenta at a 
given point in the spiral arm. As they collide though, the 
gas particl es (or small c l umps of gas) exchange angular mo- 
mentunjfl (jPobbs et al.l l2006h . The angular momenta and 
therefore velocities of the gas particles converge and locally 
a clump forms. Gas outside this locus, which has more dis- 
parate orbits, does not come sufficiently close to the density 
perturbation to experience much change in momentum. The 
locus of orbits over which gas can accumulate into a clump 
will depend primarily on the strength of the shock, i.e. the 
degree to which the orbits are perturbed by the shock. 

One measure of the strength of the spiral shock is the 
epicyclic radius. This is related to the distance and time gas 
travels along the spiral shock (together with the pitch angle 
of the spiral potential and gas pressure), which in turn de- 
termines the size of the region from which the spiral shock 
can gather material into a single clump. In Table 1, we show 
the epicyclic radius calculated at a radius of 7.5 kpc from 
the orbits of particles subject to the potential alone. Fig. 6 
shows the spacing from Table 1 plotted against the epicyclic 
radius. We see that the spur spacing tends to be about twice 
the epicyclic radius (also equivalent to the maximum ra- 
dial extent of the gas as it passes through its orbit). For 
the stronger shock, the orbits are more perturbed, so the 
epicyclic radius and spacing between the clumps becomes 
larger. However with greater thermal and/or magnetic pres- 
sure (models E and C), the spiral shock becomes weaker. 
With a weaker shock, there is less orbit crowding, and fewer 
interactions between gas particles or clumps. Interactions 

^ Particles entering the shock are at the greatest extent of their 
orbits and have lower velocities compared with gas already in 
the shock. Thus these particles tend to gain angular momentum 
whilst material already in the shock loses angular momentum. 
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Figure 3. A section of spiral arm is shown for models B and G, 
with (top) and without (lower) self gravity. The main features 
along the spiral are labelled. These are similar in each case, al- 
though with self gravity, the clumps and spurs are more coherent. 
Without self gravity these features are more dispersed and frag- 
mented. 

also produce a smaller transfer of angular momentum, since 
the difference in velocities of gas within the spiral shock, 
and compared to gas entering the spiral shock, is less. Hence 
there is a smaller locus over which particles are subject to a 
sufficient change in angular momentum to be concentrated 
into a clump, and the spacing between the clumps is smaller. 

Table 1 also includes the mass of the largest clump along 
a section of spiral arm for each model (determined using a 
clumpfinding algorithm, see Section 3.3). The typical cloud 
masses are ^ 10^ M©, which is much larger than the Jeans 
mass for the cold gas 5000 Mq). For cloud formation by 
gravitational instabilities, the Jeans analysis predicts that 
discs with larger surface densities should contain lower mass 
clouds, whilst a stronger shock will produce lower or the 
same mass clouds depending on whether the average or spi- 
ral arm surface density is considered. For agglomeration, 
both an increased surface density or a stronger spiral po- 
tential would be expected to lead to higher mass clouds, as 
found for these results (runs D, F and H). 

3.1.3 Gravitational stability of the cold gas 

As mentioned in the previous section, Q < 1 for the cold 
gas, hence the cold component is expected to be globally 
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Figure 4. This figure shows the gas which constitutes a clump 
at an earher time. The bottom panels show a clump from model 
B after 225 Myr, which contains a mass of 1.3 x 10^ M©. The 
gas which constitutes the clump is shown on the left and the sur- 
rounding material on the right. The top left panel shows the gas 
from the clump after 170 Myr (left), and again the surrounding 
region on the right. Thus the clump appears to form by the ag- 
glomeration of smaller clumps. This is in contrast to formation 
by gravitational instabilities in warm gas, where a perturbation 
develops from uniformly distributed gas. 



unstable to gravitational instabilities. Previous simulations 
indicate th at an unm agnetized disc can be unstable even 
if Q=2 (Sh ettv Qs triker 2006). Although self gravity in- 
fluences the density and morphology of clumps in the disc, 
runaway collapse does not occur in the low density mod- 
els. Collapse may be prevented in our models by supersonic 
motions, magnetic fields or insufficient numerical resolution. 
However we stress again that the size scale of gravitational 
instabilities in the cold gas in our simulations is of order 10s 
parsecs, and thus mu ch less than the scale of the features 
found in warm gas in IShettv Ostrikerl (|2006l l. Therefore, 
ignoring stellar feedback, self gravity is not likely to signifi- 
cantly change the large scale distribution of clouds and spurs 
in our models even if collapse occurs. 

To investigate the disc stability explicitly, we performed 
further simulations without the spiral perturbation. With- 
out magnetic fields, the cold component of a two fluid disc 
with E = 4 M© pc~^ is unstable. As would be expected, the 
initial velocity dispersion is not sustained on small scales 
and collapse occurs. However with a magnetic fleld of the 
same strength as model B, the disc is stable. The overall 
outcome, i.e. that runaway collapse does not occur when 
magnetic flelds are present, is therefore the same whether or 
not the spiral potential is applied. 

When magnetic flelds are included, runaway collapse is 
prevented by the overall pressure of the gas, which com- 
prises both a thermal and magnetic component. Although 
P > 1 globally, locally the magnetic pressure can be much 
higher than the thermal pressure, and since the magnetic 
flelds are fairly disordered in the clumps, the magnetic pres- 



0.8 



0.6 



o 



0.4 - 



0.2 - 





F=4% 


F=8%/ 


/ \f=i6% ; 




'l ' 1 















0.0 



0.5 



1.0 1.5 
Period (kpc) 



2.0 



2.5 



Figure 5. The Fourier transform is shown for the distribution of 
particles along a spiral arm for models B, D and F. In each case, 
the Fourier transform for the arm with the strongest signal is 
chosen (hence for the 8% and 16% cases the spacing is the largest 
of the 4 arms, since the efficiency of the shock in concentrating 
gas into large clumps will be highest to give the clearest signal). 
The period, i.e. distance between features, shifts to higher values 
as the strength of the shock increases. 
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Figure 6. The average spacing between clouds (Table 1) is plot- 
ted against the epicyclic radius (te) associated with the potential. 
The spacing increases with the epicyclic radius (measuring the 
strength of the shock), and corresponds to approximately twice 
the epicyclic radius. 



sure is relatively iso tropic. We previously found B oc 
(iDobbs Price I2OO8 ), so the magnetic fleld is strong (^100 
/xG) in dense regions. Furthermore f3 is then < 1, typically 
between 0.001 and 0.1. To illustrate the role of magnetic 
pressure, we consider the Jeans length of a local spherical 
clump. For dense gas in our simulation, with p = 10~^^ g 



cm , the Jeans length Aj = {cgTi/Gp) 



is 



0.5 pc in the 



cold gas. For the magnetic Jeans length however, the sound 
speed is replaced by the Alfven speed, i. e. Xb = {vaTr/GpY^'^ 
where VA—\f^Csl\f^ (Chandrasekharl Il96lh . This depen- 
dence on /3 increases the local Jeans length of most gas at 
this density to around 10-100 pc. Thus collapse does not 
continue below these scales, and is prevented by magnetic 
and thermal pressure. However for the 20 Mq pc~^ cal- 
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culation (model I), the densities become sufficiently large 
to overcome the pressure, hence runaway collapse halts the 
calculations. We could use a similar argument by taking the 
local value of again m odified to include magnetic fields. 
IShettv Qstrikej hood ) similarly find that collapse occurs 
earlier (within an orbit) without magnetic fields, but this is 
rather a consequence of a strong ordered field which prevents 
collapse along the magnetic field lines. 

Nevertheless, even in the lower surface density mod- 
els there is a distribution of f3 for a given density. At higher 
resolution, there may be clumps with /3 < 1 which could col- 
lapse, but are not resolved in our simulations. We performed 
simulations without magnetic fields at different resolutions. 
The peak density in the disc increases with resolution, and 
consequently at lower resolution simulations are able to run 
for slightly longer before collapse occurs. Thus the combined 
thermal and magnetic pressures primarily prevent collapse, 
but numerical resolution is also important. 

For the low surface density models, self gravity has 
only a small effect on the gas before it enters a spiral 
shock (Appendix A). Although the density increases in the 
shock, this is counteracted by an increase in the veloc- 
ity dispersion of the ^as due to collisions b etween clumps 
(|Bonnell et al.1 l2006l : IPobbs Bonnelll l2007l ) . In the most 
dense clumps though, self gravity becomes more important 
and gas reaches even higher densities. The dense clumps are 
those which are most bound, particularly as they move away 
from the spiral shock and the velocity dispersion decays. 
Clumps which become gravitationally bound during or after 
the spiral shock are then supported by magnetic and thermal 
pressure, and as mentioned above have field strengths in ex- 
cess of 100 /iG. However, as mentioned previously, gas in the 
bound clumps may undergo runaway collapse at higher res- 
olution. Furthermore, we also neglect thermodynamic pro- 
cesses and non-ideal MHD effects such as ambipolar diffu- 
sion, which are likely to promote collapse. 

As we show in Section 3.3.2, the largest GMCs which 
form in these low surface density calculations tend to be 
unbound associations of smaller clumps. The overall pic- 
ture we present is that self gravity is likely to be important 
with regards to small, dense clumps within a GMC, but it is 
the action of the spiral shock which focuses the clumps into 
a GMC (GMA (giant molecular association) may well be 
more appropriate). Although the most dense clumps within 
the GMC may become bound, supported by magnetic fields 
and supersonic velocities, the whole GMC (or GMA) is not 
necessarily bound. 

3.2 GMC formation by agglomeration and 
gravitational instabilities (high surface 
density models) 

We now consider discs with a higher surface density, where 
self gravity is likely to have a more important role in the 
formation of GMCs and their resulting properties. Fig. 7 
shows the disc structure with a much higher surface den- 
sity, of 20 M0 pc~^. The top panel shows a single phase 
simulation containing only 10^ K gas (model M), after 265 
Myr. The clumps and spurs are purely due to self grav- 
ity, as simulations of warm g as do not show such fe atures 
in the absence of self gravity (jPobbs Bonnellll2006h . The 
spacing is roughly in agreement with the predictions of 



Elmegreen k ElmegreenI (|l983l ). and not too dissimilar to 
Shettv &: Ostrikerl (|2006h . As a further test, a simulation 
with twice the disc mass was performed (model N), and the 
spacing between clumps became clearly smaller, as would be 
expected from Eqn 8. Thus o verall these simula t ions w ith 
warm gas agree with those of IShettv Ostrikerl (|2006[ ), in 
finding that GMCs form by self gravity. 

For the 20 M© pc~^ density calculations which include 
cold gas, gravitational collapse halts the calculations after 
130 Myr (model I) and 180 Myr (model L). The middle 
panel of Fig. 7 shows a calculation with cold and warm gas 
after 130 Myr (model I). The structure of the disc is differ- 
ent to single phase calculations with cold or warm gas after 
this stage. Unlike the lower surface density calculations, the 
warm gas experiences density perturbations along the spi- 
ral arms. Large coherent complexes are located along the 
spiral arms, which consist of envelopes of warm gas, and 
dense cores of cold gas. These large complexes have masses 
of 10^ M0 or more. 

The massive GMCs in the two fluid calculation, model 
I, form by a combination of self gravity, and agglomeration 
between clumps of cold gas. The structure of the disc for 
model I deviates from the low surface models B and G, with 
the presence of larger clouds along the spiral arm, whilst the 
spiral shock is less evident. The spacing between the most 
massive complexes is however compatible with either grav- 
itational instabilities in the warm gas or the agglomeration 
scenario. The complexes are separated by 1.2 kpc, which 
is similar to both Xmax for the warm gas, and the separation 
in the low surface density models, B, G and H. There is sim- 
ilarity in the structure for models I and L, so it is plausible 
that gravitational perturbations in the warm gas are seeded 
by the most dense cold clumps. At later stages of the simu- 
lation, cold clumps tend to be confined to the minima of the 
gravitational potential due to the warm gas, increasing in 
mass by the agglomeration and accretion of smaller clumps 
of cold gas. As described in the next section, the properties 
of the massive clouds in this simulation suggest that self 
gravity has a much more dominant role in determining their 
structure than for the low surface density calculations. 

With only cold gas (model L, Fig. 7 lower panel), there 
is more substructure in the gas, and gravitational collapse 
occurs on very small scales since Xmax 10-20 pc. However, 
self gravity also has a much greater role in organising the 
large scale structure than for the low surface density mod- 
els, as can be seen by comparing Fig. 7 (lower) and Fig. Id), 
which also has only cold gas, but a surface density of 4 M© 
pc~^. Thus large scale structure is present along the spiral 
arms, similar to the two fluid model (I), and again self grav- 
ity and agglomeration contribute to the formation of GMCs. 
However for the two fluid model (I), the warm gas appears 
to help concentrate the cold gas into singular large clumps, 
which smoothes out some of the structure in the disc. 

3.3 Properties of clumps formed in these 
simulations 

In this section we compare the properties of the clumps 
formed in the low and high surface density cases. These 
could potentially be compared to observations and give an 
indication of the main processes contributing to GMC for- 
mation in a particular galaxy. In particular, we consider the 
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Figure 7. The column density of the disc is shown for models M 
(top), I (middle) and L (lower) after 265, 130 and 130 Myr respec- 
tively. The gas is either all warm (top panel), half warm and half 
cold (middle panel) or all cold (lower panel). In these simulations, 
the disc mass is 5 times larger than in Fig. 1. Consequently the 
warm gas is unstable to gravitational instabilities. With no cold 
gas, the substructure in the top panel is solely due to self gravity. 
In the middle panel, both warm and cold gas are included, and 
the large complexes are presumed to form through gravitational 
instabilities in the warm gas, and agglomeration between cold 
clumps. In the lower panel, all the gas is cold, and although the 
disc appears similar to the middle panel, the coherent complexes 
present in model I along the spiral arms are largely absent, and 
there is more substructure. 



angular momentum of clouds in our models and compare 
with observations of M33. 

To select clouds, we applied a clump-finding algorithm 
to these simulations. The algorithm selects regions with sur- 
face densities of 5X1 , where E is the initial average sur- 
face density, with a resolution of 10 pc. Only clumps with 
^ 30 particles are retained (equivalent to 7500, 15000 and 
37500 M0 for the 4, 8 and 20 Mq pc~^ surface density calcu- 
lations respectively). Although we do not include molecular 
hydrogen formation in these calculations, the densest regions 
are likely to correspond to GMCs. 

3.3.1 Angular momentum of GMCs and implications for 
the ISM 

iRosolowskv et al] (|2003 ) have determined the angular mo- 
mentum from clouds in M33 and compared these values 
with those expected from several theories of GMC forma- 
tion. Their values of angular momentum are lower than ex- 
pected if GMCs form by gravitational instabilities. Further- 
more 40% of their clouds show retrograde motion, which can- 
not be explained if self gravity is predominant. They suggest 
that magnetic braking may play a role in removing angular 
momentum from GMCs as they form. 

Fig. 8 shows the specific angular momentum of the 
clouds found in models B and I (with respect to the cloud's 
centre of mass) . We see that in both cases the GMCs exhibit 
prograde and retrograde rotation. The distribution of angu- 
lar momentum and ratio of prograde and retrograde clouds is 
similar in models B and G, i.e. with and without self gravity, 
with approximately twice as many prograde clouds as ret- 
rograde. Retrograde clouds arise through collisions between 
clumps (see also Tasker et al. 2008), and the frequency of 
retrograde clouds which arises in these models is an indi- 
cation of the clumpy nature of the gas in the simulations. 
Retrograde motions are not expected if the gas is treated 
solely as a fluid, and accordingly all the GMCs (formed by 
gravitational instabilities) in model M (with only warm gas) 
were found to exhibit prograde rotation, although there are 
far fewer objects. Thus observations of retrograde clouds in 
M33 suggest that at least some component of the ISM is 
clumpy. The magnitude of the specific angular momentum 
in models B and G, where clumps form mainly by agglom- 
eration, does not exceed 100 km s~^ pc an d is co mparable 
to the values obtained by 'Rosolowsky et alj (|2003l l for M33. 

When the surface density is 20 Mq pc~^, the distribu- 
tion of clouds of mass ^ 4 x 10^ Mq is similar, and again 
there are twice as many prograde clouds as retrograde. How- 
ever there is a surplus of prograde clouds with a high an- 
gular momentum. These are the most massive clouds where 
self gravity is significantly contributing to the growth of the 
cloud. The rotation of the clouds is visibly evident from the 
presence of spiral tails of gas (Figs. 7 and 10). Again these 
most massive clouds are associated with instabilities in the 
warm component of the ISM, which behaves as as fluid. 

3.3.2 Mass spectrum and gravitational stability of the 
GMCs 

We also compared the boundedness of clouds in these calcu- 
lations. As expected, in the simulations with a higher surface 
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density, the clumps are more bound. For model B, where 
S = 4 M0 pc~^, the virial ratio (a) is between 1 and 6 for 
the clumps of mass > 10^ M©, so these clouds are unbound 
{a — ba^ R/3GM where a is the 3D velocity dispersion, R 
is the radius of the clump and M the mass) . Any marginally 
bound clumps in model B are sufficiently supported by mag- 
netic and thermal pressure for the calculation to continue. 
With model I, where E = 20 M© pc~^, a lies between 0.7 
and 2 for clouds of mass > 10^ M©. In model I, the magnetic 
and thermal pressure are not sufficient to prevent collapse 
in the most bound clouds, and the calculation stops. In the 
low surface density runs, the clouds are more bound when 
self gravity is included. Without self gravity (model G), a 
is distributed about a mean of approximately 5 for clumps 
> 10^ Mo. 

The distribution of clouds also changes with self gravity, 
as the simulations with self gravity contain more interarm 
clouds. Clouds formed in the spiral arms are more bound, 
and remain intact for longer in the interarm passages. We 
also found that the clouds in model I, with the higher sur- 
face density, produced a shallower mass spectrum (Fig. 9), 
with dN/dM oc M~^-^^ compared to for model B (self 

gravity and low surface density) and M~^'^ for model G (no 
self gravity). At high masses, self gravity increases the ag- 
glomeration of clumps and accretion onto clouds, so there 
are more massive clouds and a shallower slope. A similar re- 
sult was produce d recently in simulations of prestellar cores 
(|Dib et al.ll2°Q0i ). 

Fig 10 shows individual clouds from the low and high 
surface density calculations, B and I. As mentioned in Sec- 
tion 3.1.2, clouds formed by agglomeration also tend to con- 
sist of numerous clumps dispersed throughout the GMC. By 
contrast, in the high surface density run where the GMCs 
are more bound, the GMCs are much more centrally con- 
centrated. 

Generally, the properties of clouds formed in model I, 
where self gravity plays a strong role in the formation of 
GMCs are most similar to those observed for the Milky Way 
(i.e. lower a and the shallower mass spectrum). Even the ve- 
locity sizescale relation of the clouds resembles the a oc r^'^ 
observed law, whilst there is no clear relation in the other 
simulations. However, the properties of clouds in model I 
are dissimilar to M33, where there are no 10^ Mp G MCs. 
The mass spectrum is also steeper telitz e t "aD 12007 ^) . and 
the GMCs all have relatively low specific angular momenta. 



4 DISCUSSION 

The formation of molecular clouds in spiral arms is be- 
lieved to occur primarily by gravitational instabilities or 
cloud collisions. For two decades however, agglomeration 
models for molecular cloud formation have been largely ne- 
glected. We have improved previous computational analysis 
of cloud format i on via collisions (e.g. Casoli &: Combesll 19821: 
TomisakalfTisi : iHausman k Robertslll984l : iKwan Valdes" 
19871 ) by being able to perform full hydrodynamical simu- 
lations, which do not require clump masses or the degree 
of fragmentation or coalescence as input parameters. Fur- 
thermore, in our models, the constituent clumps which form 
the GMCs are generally not self gravitating, avoiding the 
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Figure 8. The specific angular momentum is plotted versus the 
mass of a cloud for the low surface density case (model B), top, 
and the high surface density case (model I). Self gravity is in- 
cluded for both. With the higher surface density though, the dis- 
tribution of angular momentum of the clouds extends to over 100 
pc km s"-*^. These high angular momentum clouds are formed 
mainly by self gravity and exhibit prograde rotation. 
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Figure 9. The mass spectrum is shown for clouds extracted from 
Runs G (red), B (green) and I (blue). The mass spectrum becomes 
shallower as self gravity becomes more significant. The shallowest 
slope is obtained when the surface density is 20 pc~^. 
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Figure 10. These figures zoom in on individual structures in the disc which would correspond to GMCs. Vectors show the magnetic 
field integrated through the plane of the disc. The left panel shows a cloud from the lower surface density model (B), where E = 4 
pc"-^. The mass of the cloud is 3 x 10^ M©, which is distributed in several separate clumps. The cloud on the right panel contains ~ 10^ 
M0 and is taken from model I, with the 20 Mq pc~^ surface density, and self gravity has much more influence on the structure of the 
cloud. The rotation of the cloud has produced tidal tails extending from the GMC, and the mass of the cloud is concentrated in the 
centre of the cloud. 



timescale problems associated with previous models. How- 
ever, magnetic pressure is also important in the absence of 
feedback, especially as more massive clouds develop in the 
spiral arms. 

In the lower surface density results we have presented, 
the formation of GMCs is due to orbit crowding in the spiral 
shock, which leads to the agglomeration of clumpy gas into 
larger structures. Without the spiral shock, interactions be- 
tween clumps are much less fr equent and collision s predomi- 
nantly lead to fragmentation toittins et ahlboosl ) . Although 
fragmentation occurs in our models, the clumps in the spi- 
ral arms tend to be more massive and dense than the mate- 
rial entering the shock, so fragmentatio n is reduced. Orbit 
crowding was previously suggested by [Roberts StewartI 
('1987') as the dominant means of GMC formation in spi- 
ral galaxies. Their calculations allowed dissipative collisions 
between clouds, rather than assuming that colliding clouds 
coalesce into coherent objects, as was the case for other mod- 
els at the time. This condition allowed structure to emerge 
along the spiral arms, similar to our hydrodynamical mod- 
els. We also find that GMCs formed by orbit crowding tend 
to be associations of smaller clouds rather than gravitation- 
ally bound structures, again in agreement with arguments 
presented in [Roberts StewartI ([l987). 

Stability analysis of the ISM subject to gravitational 
collapse has been successful in explaining the observed sep- 
aration and masses of GMCs along spiral arms. No such 
analogous interpretation has previously been identified for 
collisional models. The Jeans analysis however does not ex- 
plain the spacing between spiral arm clouds in our lower 
surface density simulations. Instead we interpret this spac- 
ing in terms of orbits of the gas and relate the distance be- 
tween GMCs to the epicyclic radius. We obtain separations 
of approximately 1-1.6 kpc (although dependent on the pa- 
rameters for the spira l potential), which is at the lowe r end 
of those observed by lElmegreen ElmegreenI (|l983h , but 



iLa Vigne eraP (|2006l l find spacing between feathers some- 
what smaller than this (200 pc -1 kpc depending on galactic 
radius). For the agglomeration scenario, the mass of GMCs 
increases for stronger shocks or a higher surface density. 
Although observations of molecular cloud masses are lim- 
ited in external galaxies, M 51 contains GMCs of > 10^ 
M© (iRand Kulkarnilll990l ) whilst galaxies without strong 
spiral shoc ks, such as M31 and M33 contain no GMCs of 
> 10^ Mo (Sheth et al. 2000VEngargiola et al." 2003). How- 
ever, more massive clouds tend to be formed anyway when 
the disc is gravitationally unstable. 

There are several possible issues with the physics as- 
sumed for these calculations. Although we state the value 
of (3 in these simulations, in the calculations with cold gas, 
the magneti c field tends to be low er than observations, as 
described in " Pobbs Pried (|2008l l . This is because we as- 
sume an initially uniform density distribution, whereas the 
cold gas in the ISM is situated in dense clumps. However, 
this is unlikely to effect our conclusions, since the ISM is ob- 
served to be clumpy despite higher field strengths, and this 
is the main requirement for agglomeration to occur. Nev- 
ertheless, calculations in which the magnetic field strength 
is consistent with observations in both the warm and cold 
component are needed in future work. 

An important caveat with our results is that we have not 
included heating and cooling of the ISM. Our models require 
that cold gas exists prior to entering the spiral shock. Obser- 
vations show a substantial part of t he ISM in the Milky Wa y 
is located in cold dense HI regions (jHeiles TrolandH2003l ) . 
though it is unclear whether this is ubiquitous through the 
disc, or just associated with gas that is already in the process 
of forming molecular clouds. If it is the latter, then GMC 
formation may occur primarily by thermal and gravitational 
instabilities, with gas cooling and regions fragmenting as col- 
lapse occurs. By including a detailed model of the thermo- 
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dynamic processes of the ISM in our upcoming calculations, 
we will be able to address this issue. 

Furthermore, we have used the thermal sound speed of 
the gas, rather than the r.m.s. speed associated with turbu- 
lence in the gas for the Jeans analysis of the cold gas. From 
movies of gas flowing through the spiral arm, it is evident 
that the gas flow is fairly chaotic, and thus the velocity dis- 
persion is somewhat higher than the thermal sound speed 
of the cold gas. Nonetheless we do not get kpc size clumps 
supported by turbulence. Our previous simulations of star 
forming cores also showed that turbulence does not act as an 
isotropic pressure supporti ng the cloud against gravitational 
collapse (|Dobbs et al.l 1200 5). Lastly, we do not include stel- 
lar fe edback, which will disrupt structure along the spiral 
arms (|Wadall2008l : IShettv Ostrikei]l2008l ). 

In addition to this further physics, direct comparisons 
with actual galaxies will be required to evaluate which are 
the dominant mechanisms for the formation of GMCs. Inter- 
estingly, the surface dens ity of our Galaxy is around 10 
within a radius of 12 kpc (|Wolfire et al.ll200"3l ), roughly at the 
point gravitational instabilities start to influence the large 
scale structure in our models. 



5 CONCLUSION 

We have investigated GMC formation by agglomeration and 
self gravity. Agglomeration occurs in spiral shocks providing 
there is a clumpy constituent of the ISM, assumed to be 
cold HI or molecular gas. This process occurs regardless of 
whether self gravity is included, and in low surface density 
calculations, the formation of GMCs is predominantly due to 
agglomeration. The converse situation arises when there is 
only warm gas, and GMC formation occurs by gravitational 
instabilities, providing the disc is gravitationally unstable. 
Generally, both agglomeration and self gravity are expected 
to contribute to GMC formation, with self gravity becoming 
more important to GMC formation and disc structure as the 
surface density increases. 

The degree to which these processes determine GMC 
properties will depend on the surface density of the galaxy, 
the thermal nature of the ISM, and most likely the mag- 
netic field strength. In particular, GMCs with retrograde 
rotation can be produced when the ISM is clumpy. A caveat 
to these calculations is that we assume a two fluid model 
with no heating or cooling. We expect to include the ther- 
modynamics of the ISM in future work. A further caveat 
is that gravitational collapse should occur in our models 
and lead to stellar feedback. Cooling of the gas to temper- 
atures < 100 K, the inclusion of non-ideal MHD processes, 
and possibly higher numerical resolution would induce col- 
lapse. However the aim of the current paper is to investigate 
the structure of the disc without feedback. The aim of fu- 
ture simulations will be to see how this picture changes with 
stellar feedback. 
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APPENDIX A: CLUMPINESS OF THE GAS 
PRIOR TO A SHOCK 

In order for agglomeration to occur in the shock, the ISM is 
assumed to be clumpy. The ISM is observed to be inhomo- 
geneous over many scales, and C NM (cold neutral med ium) 
features are observed in 100 pc (IStrasser et al ]l2007|) and 
parsec/subparsec structures fe.g. lPickev Lockmanlll990l V 
The dumpiness of the gas prior to the shock in our mod- 
els arises through the Poisson noise of the particle distri- 
bution. We show the density distribution of particles be- 
fore gas has entered the spiral shock in Fig. Al (top), for 
model B. The particles are selected from a 1.5 by 1.5 kpc 
region after 40 Myr. The particles are about to enter a spi- 
ral shock, but at this early stage in the simulation, have not 
been significantly perturbed by any previous shock. At the 
very beginning (t=0) there is only an order of magnitude 
variation in density. However the gas is then subject to the 
disc potential (Eqn (1)), hence the density varies with scale 
height with low density gas at high z. Fig. Al also shows the 
smoothing length, h, of the particles. So we see that the den- 
sity spans nearly 5 orders of magnitude, with corresponding 
length scales of a few to >100 pc. 

The density distribution is also shown for different mod- 
els in Fig. Al (lower). The density distribution is narrower 
for the single phase case, without the warm, lower density, 
component. With self gravity, the densities are only slightly 
higher indicating self gravity has little effect on the distri- 
bution (self gravity only effects a small fraction of the mass 
which occupies high densities) before the gas enters the spi- 
ral shock. Overall any differences in the distributions for 
models B, G and J have little impact on the large scale struc- 
ture that emerges because the size scales of the initial struc- 
ture is always much less than the sizes of the large clouds 
that emerge later in the simulations. Finally in Fig. Al, the 
density distribution is shown for an interarm region at a 
much later time in model B. The much higher densities in- 
dicate the influence of self gravity in the spiral shock. The 
density of the gas increases in the shock such that self gravity 
has much more of an effect. Hence clumps with much higher 
internal densities occur at later times in the simulation (and 
compared to the non-self gravitating model). 
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Figure Bl. This figure shows the disc and a subsection of the disc from model J, with 1 miUion particles (left) and 8 million particles 
(right). The equivalent simulation with 4 million particles is shown in Fig. Id). The lower panels intend to show that the location of large 
scale features is similar in each simulation and in particular, the number of large clumps per unit length along the spiral arm does not 
change. However there is clearly much more substructure in the higher resolution calculations and interarm spurs are less well defined. 
The inner arm appea rs more continuous because the spiral potential is weaker at these radii (the strength of the spiral potential from 
ICox Sz Gomed (|2002l ) has a maximum at R=8 kpc). 



APPENDIX B: RESOLUTION STUDY 

We investigate the difference in structure depending on res- 
olution by repeating model J with 1 and 8 million particles. 
The results are shown in Fig. Bl, which displays the whole, 
and a subsection of the disc. The structure of the disc clearly 
changes with resolution. The interarm spurs are less distinct 
and more disjointed in the 8 compared to 4 million parti- 
cle simulation. However the actual location and number of 
features along the spiral arm appear similar at different res- 
olutions. In particular the number of large clumps per unit 
length is unchanged. Although some of the spiral arm clumps 
are much easier to distinguish than others: structures such 
as that labelled '6' are instead collections of many smaller 
clumps. The similarity of spacing between the clumps is ex- 
pected, since the size scale of the initial dumpiness in the 



gas is much smaller than the resulting size scale of clouds 
and spurs which form in the disc. However for higher reso- 
lution, it is more difficult to assemble gas into a coherent, 
single cloud. 



Overall, the increased resolutio n has a similar e f fect t o 
a previous resolution comparison in iDobbs k Prk3 (|2008h . 
in that the same large scale features occur in different res- 
olution simulations, but the structure is more disjointed 
at higher resolution. In the calculations here, the gas also 
reaches higher densities in the 8 million particle case (com- 
pared to 1 million), since the dense regions of GMCs are 
better resolved. 
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Figure Al. These panels consider the density distribution of 
the gas just before entering a spiral shock. Particles are selected 
from a 1.5 x 1.5 kpc interarm region after 40 and 265 Myr during 
models B, G and J. The top panel shows the number of particles 
at different densities (crosses) in model B, the two fluid case with 
self gravity. The smoothing length (stars) gives an indication of 
the length scale associated with the structures. The lower panel 
shows the distribution of particles in the same region for models 
B (two fluid case with self gravity, triangles), G (two fluid case 
but no self gravity, stars) and J (100 K gas only, diamonds). 
There is little difference with or without self gravity. Finally the 
density distribution is shown at a later time of 265 Myr for model 
B (squares). The gas has reached much higher densities at this 
point due to self gravity acting on dense gas in the spiral shock. 
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